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Abstract. On the basis of generalization of upper and lower solution method to the 
singular two point boundary value problems, the existence theorem of solutions for the 
system, which models a process of magnetic insulation in plasma is proved. 
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1. Introduction. 

Investigation of mathematic models of magnetic insulation has been started by P.Degond, 



N.Ben Abdallah and F.Mehats in 1995 year. In 1996 P.Degond has put to the author of this 



paper the problem on existence of solutions of limit system (I) and its generalization to the 



problem with free boundary. The effect of magnetic insulation consists in that the electrons 



emitted from cathode cannot reach the anode due to the extremely high applied electric and 
magnetic field; they are reflected by the magnetic forces back to the cathode. Thus there 



is electronic layer outside of which electromagnetic field is equal to zero [1] . Here two basic 
regimes are possible: the first, when electrons reach the anode — " noninsulated" diode and 
the second one, when electrons rotate back to the cathode — "insulated" diode. The regime 



of "noninsulated" diode is described by the following nonlinear two-point boundary value 
^ ■ problem 



d 2 LD 1 + (p(x) A 

d% 2 ^ (I + p{x)) 2 - I - a{x) 2 



(I) 



„ - j x — j= = == = G((p, a); a(0) = 0, a(l) = a L , 



d 2 a . a(x) 

d % 2 X ^J(l + ip(x)) 2 -l-a(xy 

where j x > 0, x G [0, 1]; <p is the potential of electric field and the potential of magnetic field 



is a. 



Our main goal consists in search of positive solutions of system (I) that is if > 0, a > 
and their dependences upon parameter j x . Here there are some interesting questions about 
solvability of this problem, because the system (I) is singular in zero for ip = and in this 
connection, we can not say about properties of monotonicity of right parts on the interval 
(p G [0, oo) and, hence, about Lipschitz condition. The problem (I) has no a property of 
quasimonotonicity in cone. Thus a standard upper and lower solution method, developed 
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for the systems of semilinear elliptic equations in partially ordered Banach space [3], does 
not work. In spite of this fact, we show the existence of lower and upper solutions of problem 
(I) without conditions of local Lipschitz continuity and quasimonotonicity using sufficiently 
simple techniques. To this purpose, we modify the P.J.McKenna and W.Walter [4] theorem 
of existence of lower and upper solutions for arbitrary elliptic systems 

Au + f(x,u) = f2, u = on dVl. 

where u = (ui, . . . , u n ), f = (/i, . . . , f n ) are n — vectors, Q — is an open bounded subset of 
R M with smooth boundary dQ, and f(x,u) is uniformly Holder continuous (with exponent 
a) in x and Lipschitz continuous in u. 

In section 2 we will prove theorem 2 and propositions 1, 2 on the existence of semitrivial 
solutions of problem (I) by upper and lower solution method. The estimations to the value 
of electrostatic potential on the anode y?L and the current j x are obtained. In section 3 we 
formulate the principal theorem 4 on the existence of positive solutions of problem (I) and 
the estimation to the value of magnetic field on the anode <zl is given. 

We note that system (I) was studied in [2] by a shooting method with (3 = a'(0) and j x 
as shooting parameters. The strategy is: given the values of (3 and j x , solve (I) with the 
Cauchy conditions ip(0) = 0, a(0) = 0, <p'(0) = 0, o'(0) = f3, and then adjust the values in 
order to fulfill the conditions ip(l) = tp L and o(l) = aj_,. 

2. Existence of semitrivial solutions of problem (I). 

Let us introduce the definition of cone in a Banach space X. 

Definition 1. Let X be a Banach space. A nonempty convex closed set P C X is called 
a cone, if it satisfies the conditions: 
(i) x G P, A > implies \x G P; 

(ii) x G P, —x G P implies x — O, where O denotes zero element of X . 
< is the order in X induced by P, i.e. x < y if and only if y — x is an element of P. 

We will denote [x, y] the closed order interval between x and y, i.e. 

[x,y] = {zeX :x<z<y}. (1) 

We will also assume that the cone P is normal in X, i.e. order intervals are norm bounded. 

InX 

X = {(u,v) : u,v G C 1 (Cl),u = v = 0} 

we introduce the norm \U\x = M^i + \v\c^, and the norm \U\x = Moo + Moo i n C, where 
U = (u,v). Here P is given by 

P = {(u,v) e X : u > 0, v > for all x G Q}. (2) 

So, if u 7^ 0, v t^ belong to P, then —u, —v does not belong. We will work with classical 
spaces on the intervals I — [a,b], I —]a,b], I = (a, b): 



C(I) with norm || u ||oo= max{\u(x)\ : x G /}; 
ai(T\ —II 11 II -i- II 11' II • 

\ A ) ~\\ U lloo T || " ||oO! 

Ci oc (I), which contains all functions that are locally absolutely continuous in /. We introduce 
a space Ci oc (I) because the problem (I) is singular for (p = 0. The order < in cone P is 
understood in the weak sense, i.e. y is increasing if a < b implies y(a) < y(b) and y is 
decreasing if a < b implies y(a) > y(b). 

Theorem l.[5] (comparison principle in cone) Let y G C(I) f]Ci oc (I). The function f is 
defined on I x R. Let f(x,y) is increasing in y function, then 

v" — f(x, v ) > w" — f(x, w) a.e. on I, (3) 

v(a) < w(a), v{b) < w{b) 

implies 

v < w on I. 

Remark 1. Let f(x,y) is decreasing, then theorem 1 remains without changes, if both parts 
of (3) multiply onto -1. 

For the convenience of defining an ordering relation in cone P we make a transformation 
for the problem (I). Let F(ip,a) and G(f,a) be defined by (I). Then throuth the transfor- 
mation ip = —u the problem (I) is reduced to the form 

d 2 u 1 — u A f,, 

~T^ =3x I, xo == 7 = F ^x,u ) a), w(0) = 0, u(l) =<p L , 

ax w(l — u) 2 — 1 — a 2 

(II) 

— = j x = G(j x , u, a), o(0) = 0, o(l) = a L . 

ax \ (1 — u) 2 — 1 — a 2 



We note that all solutions of the problem (I), as well the problem (II), are symmetric with 
respect to the transformation of sign for the magnetic potential a : (</?, a) = (</?, —a) or 
the same (u, a) = (u, —a). Thus we must search only positive solutions (p > 0, a > in 
cone P or only negative ones: <p < 0, a < 0. Thanks to the symmetry of problem it is 
equivalently and does not yields the extension of the types of sign- defined solutions of the 
problem (I) (respect. (II)). Once more, we note that introduction of negative electrostatic 
potential in problem (II) is connected with more convenient relation between order in cone 
and positiveness of Green function for operator — u" that we use below. 

Definition 2. A pair [(<p ,ao), (ip°,a )] is called 
a) sub-super solution of the problem (I) relative to P, if the following conditions are satisfied 

(<p , a ) g c loc (i) n c(i) x c loc (i) n c(i), 

- ; (4) 

(y?°, a ) G C loc (I) C{I) x C loc (I) f] C(I) 
V? - j x = F((p , a) < in /, 



(5) 



„ov/ l + V 



iv) ~ jx > = FUp", a) > in / Va e [a , a u 

^1 + y? ) 2 - 1 - a 2 



Oq - j x . = G(<f, a ) < in J, 

^(l + ^-l-a 2 



( a °)"- Ja — = a ° = G(<p,a ) > in / V^ e KA 



(6) 



V(i+ ^ - 1 - (o°) 2 

Po < <^°, ao < a in / (7) 

and on the boundary 

y?o(o) < o < </(0), v?o(i) <^l< y°(i), 

(8) 
a (0) < < o°(0), a (l) < «l < a°(l); 

b) sub-sub solution of the problem (I) relative to P, if a condition (4) is satisfied and 

<Po ~ F(j x , tpo, ao) < in I, 

(9) 
Oo -G(j x ,ip , o ) < in J 
and on the boundary 

V?o(0)<0, v^o(1)<^l, a (0)<0, ao(l)<o L . (10) 

Remark 2. In definition 2 the expressions with square roots we take by modulus |(1 + (p) 2 — 
l-a 2 |. 

By analogy with (9), (10) we may introduce the definition of super-super solution in cone. 

Definition 3. The functions $(x, x ai ,j x ), ^i(x,x (pj ,j x ) we shall call a semitrivial solu- 
tions of the problem (I), if $(x, x ai , j x ) is a solution of the scalar boundary value problem 

<f" = F(jx, f, x^) = j x h = , (HI) 

1 + V?) 2 - 1 - (x a J 2 



<p(0) = 0, <p(l) = <p L , 
and ^(xjX^jjx) is a solution of the scalar boundary value problem 

a" = G(j x , x Vj , a) = j x (IV) 

,1 + x Vj ) 2 - I- a 2 



a(0) = 0, o(l) = a^. 

Here x ai , i = 1, 2, 3 and x Vj , j — 1,2 are respectively, the indicators of semitrivial solutions 

$(x,x ai ,j x ), ^i(x,x Vj ,j x ) defined by the following way: 

x ai = 0, if a(x) = 0; 

x a2 — a ° 1 if a = a° be upper solution of the problem (IV); 

Xa 3 = 0,0, if a = ao be lower solution of the problem (IV); 

x Vl — tp°, if tp — tp° be upper solution of the problem (III); 

x ip 2 = ¥0, if V 9 = ¥0 be lower solution of the problem (III). 

From definition 3 we obtain the following types of scalar boundary value problems for 
semitrivial (in sense of definition 3) solutions of (I) (resp.(II)): 

<p» = F(<p,0)=j x l + ^ , y?(0) = 0, <p(l)=<p L . (AO 

^/(1 + V^) 2 - 1 

<p" = F(<p,a°)=j x 1+ip <f(0) = 0, <p(l) = VL. (A 2 ) 



^(l + y^-l-K^ 



^=^,00)=^ ,, * + v? = = , y?(0) = 0, ¥»(1) = ¥»l. (A3) 

;i + v?) 2 -i-(a ) 2 



a" = G(tp°, a) = j x a(0) = 0, o(l) = a L . (A 4 ) 

yj(l + v? ) 2 - 1 - a 2 

a" = G(y?o, a) = j x , =, a(0) = 0, o(l) = a L . (A 5 ) 

^(l + y^-l-a 2 

We shall find the solutions of problems (A\) — (A3) with condition 

¥0 < <p°, 

where (po(x ai ), tp°(x a2 ) are respectively, lower and upper solutions of problem (Ai). The 
solution (tp, a) of problem (I) should be belong to the interval 

tp e$(tp,0)f)$(tp,a o )f)<S>(tp,a o ), 

ae $i(v? ,a)fl $ i(^o,a). 
Moreover, the ordering of lower and upper solutions of problems (A\) — (A3) is satisfied 

tp (x ai ) < tp (x a2 ) < tp (x a3 ) < tp°(x a2 ) < tp°(x ai ). 

We shall seek the solution of problems (A4) — (A 5 ) with condition 

a < a , 

in this case the following ordering of lower and upper solutions of problems (A4) — (A$) 

oo(^i) < 00(^2) < o°(x ¥ , 2 ) < o (x Vl ). 



is satisfied. 

We go over to the direct study of the problem (III), which includes a cases (A^ — (A 3 ). 
Let us consider the boundary value problem (III) with 

F(x,y9):(0,l]x (0,oo)^(0,oo). (B x ) 

In condition (Bi) for F(x, tp) we dropped index Oj, considering a general case of nonlinear 
dependence F of x. 

We shall assume that F is Caratheodory function, i.e. 

F(-,s) is measurable for all s G R, (B 2 ) 

F(x, •) is continuous a.e. for x G]0, 1], (B 3 ) 

and the following conditions hold 

/ s(l - s)Fds < 00. (B 4 ) 

JO 

OF /dip > 0, i.e. F is increasing in p. (B 5 ) 

There are 7(2;) G -^QO, 1]) and a E R, < a < 1 such that 

\F(x,s)\ < 7 (x)(l + |s|- a ), \/(x,s)E]0,l]xR. (B 6 ) 

We are interested in a positive classical solution of equation (III), i.e. </? > in P for 
x G]0, 1] and ip G C([0, l])f|C 2 (]0, 1]). The problem (III) is singular, therefore, condition 
(Bi) is not fulfilled on the interval tp G (0, 00) and in this connection, the well-known 
theorems [3] on existence of lower and upper solution in cone P does not work. It follows 
from theorem 1, since F in (III) is increasing in tp, then tp < w for x G]0, 1], where tp and w 
satisfy to the diffeential inequality (3). 

Theorem 2. Assume conditions (B2) — (Bq). Then there exists a positive solution 
tp G C([0, 1]) f|C 2 (]0, 1]) to the boundary value problem(III) . 

Proof. Let tp > is a solution of problem (III). By theorem 1 tp < w for x G]0, 1]. Take 
e > and consider equation 

V? e = i* 1 NO = == = ^e(jx, </?e + e, x a J. 

1 + ^ + e) 2 - 1 - (x ai ) 2 



(11) 

V9 e (0) = 0, ^(1) = ^L- 

Let w and </? are upper and lower solutions of equation (11) (below, in proposition 1 is shown 
that such solutions really exist). Hence the theorem on monotone iterations [6] gives an 
existence of classical solution tp e of equation (11), which satisfies w > tp e > tp for x G]0, 1] 
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and is bounded in C. Thus F e (j x ,tp e + e,x ai ) is bounded and there exists uniform limit 
lim^o^e — V 9 - ^ follows from the last, if < 77 < \, then \im e ^ F e (j x ,tp t + e, x 0i ) = 
F(j x , tp, x ai ) uniformly on [77, 1 — 77] and tp > for x G [77, 1 — 77]. 

Since y2 e converges uniformly on [0, 1], this is implies the existence lim e ^ ¥e(v)- Therefore 
there exists lim £ _+o ty? e (x) on compact subsets (0,1) and {tp e } uniformly converges on (0,1) to 
the different iable function tp on [77, 1 — 77]. It follows from the last, tp is twice differentiate 
on [77, 1 - 77], tp" = F(j x ,tp,x ai ), x G [77, 1 - 77] and u G C([0, l])nC 2 (]0, 1]) is a positive 
solution of problem (III). 

Remark 3. Delicate moment in proof of theorem 2 is connected with finding of a lower 
tp and an upper w solutions for perturbated problem (11). As a lower solution we can 
take solution of equation (A±) (semitrivial solution tp), then an upper solution will be, for 
example, maximal solution of equation (-Ai). 

Application of monotone iteration techniques to the equation (III) gives an existence of 
maximal solution tp(x,j x ) such that 

tp(x, Xj) < tp(x, Xj) < w(x) x g]0, 1]. (12) 

Proposition 1. Let < c < j x < j™ a:r . Then equation (Ai) 

1, 1 + tp 
tp = F(j x , tp, 0) = Jx ==, 

V^(2 + <p) 

<p(0) = 0, <p(l) = <p L 

has a lower positive solution 

u = (5V/ 3 , (13) 

if 

4<5 3 > 9j™ a:E (l + 5 2 )/V2TP (14) 

and an upper positive solution 

u° = a + (3x (a,(3>0), (15) 

here 

Vl > S 2 , (16) 

where 5 is defined from (14). 



Remark 4. Square root is taking as J\tp(2 + tp)\ in the case of negative solutions. Here 
u° = —ex is an upper solution, and Mo = _ 2 + e is a lower solution (0 < e < 1). Hence, 
equation (A ± ) has negative solution only for < tp L < — 2 because F(x, —2) = —00. 

It follows from (14), (16) that a value of current is limited by the value of electrostatic 
potential on the anode tp L 

j x <jr x <Hv L ). (17) 
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Analysis of lower and upper solutions (13), (15) exhibits that for 5 2 = tp^ > 2 and a — /3 < 1 
interval in x between lower and upper solutions is decreased, and for the large values of 
potential </?£, diode makes on regime ^x 4 / 3 . 

Proposition 2. Let < c < j x < j™ ax . Then equation (A4) 



G(j x ,ip°,a) =j a 



l + (p 



0A2 



I -a 2 



o (0) = 0, a(l) 



a L 



with a lower solution clq = and an upper solution a = u° > 0, conditions (14), (16) has 
an unique solution a(x,j x ,c), which is positive, moreover 



0<a L < vV (2 + tp) 



(" 



Proof. The positive solution of problem (Ai) is concave and be found as a solution of initial 
problem with a(0) = 0, a'(0) = c, where c is shooting parameter. The solution a = a(x, j x , c) 
is unique and strongly decreasing in c because right part of differential equation is decreasing 



in a. The least nonnegative solution is f(x,j x , 0) = and for < ai < J(p° L {2 + ip° L ) exists 
only one solution and no positive solutions for other values a^. 

Remark 5. The problem (A 5 ) is considered by analogy with problem (A 4 ), change of an 
upper solution a = w c 



to a lower a = u one and < a L < Jip 0L (2 + <^ol)- 

Following to the definition 2 and propositions 1,2, solutions of problems (III), (IV) we 
can write in the form (fig. 1): 



if, a 




I x 

Figure 1: location of lower (ip , a ) and upper (ip°, a ) solutions 

lower-lower (ip , a )): 

(po = u = S 2 x 4/3 , a = 0, <p L > 5 2 ; 

upper-lower (ip°,a Q ): 

ip° = u° = a + (3x, ao = 0, 5 2 < ifL < C, C = max{a, /3}; 



lower-upper (ip , a°): 

ip = u = 5 2 x^ 3 , a° = u°, (f L >S 2 , a L < ^(u (2 + u ); 
upper-upper (ip°,a°): 

p = u ° = a + /3x, a = u°, y? L <C,a L <a° < u°. 

3. Existence of solutions of system (I). 

In the previous section we demonstrated the existence of semitrivial solutions of system 
(I). Here we show existence of solutions for the complete system (I), using the following 
McKenna- Walter theorem. 

Theorem 3. [4] Assume conditions {Bi) — (Bq). We assume that there exists the ordered 
pair (u,u) — lower and upper solutions, i.e. 

u,ueQ oc ((0,l]) 2 f]C([0,l]) 2 , u<u ]0,1] 

m(0)<0<m(0), u(1) < u L < w(l); u L = { VL ,a L ), 

Vx 6]0, 1] : Vz e R 2 , 

u(x) < z < u[x\ Zk = u k (x); 

—lLk( x ) > h k (x,z) (19) 



and 



?2 



Vx e]0, 1] : Vz G R 2 
u(x) < z < u(x), Zk = Uk(x) : 

—u k (x) < h k (x,z) (20) 

for all k G {1,2}. Then there exists a solution u G C 2 ((0, l]) 2 f|C([0, l]) 2 of the problem 

- U " = h(;u(-)) ]0,1] 

u (0) = 0, «(1) =u L . 

For keeping of ordering of lower and upper solutions in theorem 3 (in cone P) we write 
differential inequalities (19), (20) in the following form 

\/z G [v(x),w(x)], z\ = Wi(x) : 

± w'[(x) < ± F 1 (w 1 (x),z 2 ) 
\/z G [v(x),w(x)], z\ — Vi(x) : 

± «i(x) I ± F 1 (v 1 (x),z 2 ) 



\/z G [v(x),w(x)]] Z 2 = w 2 (x) 

± w' 2 \x) < ± F 2 (z 1 ,w 2 ) 
Mz e [v(x),w(x)\; z 2 = v 2 (x) 

^(x) > ± F 2 (z 1 ,v 2 ). 

Remark 6. Change of signs with (+) to (— ) in differential inequalities is connected with 
adjustment of signs and ordering (<) of lower (upper) solutions of system (I) in definition 2 
and lower (upper) solutions in theorem 3. 

From the last relations we obtain 

w"(x) = Fi(wi(x),0) < F 1 (wi,z 2 ) 

(21) 
v'i(x) > supF 1 (v 1 (x),z 2 ) 

w 2 (x) < F 2 (Z!,W 2 ) 
v 2 (x) > supF 2 (z 1 ,v 2 ) 

From inequality v 2 (x) > sup 2i F 2 (zi,v 2 ) we get estimations to the value of magnetic field on 
the anode a^ 

aL< h < ir: < ^£il (23) 

2 ~ 2 ~ 2 y J 

taking account of (17) and ©^ > 0. Under realization of estimation (23) diode works in 
" noninsulated" regime, moreover the value a^ is limited by value of electrostatic potential 
on the anode </?£ with critical value </?£ = 2. In increasing of magnetic potential a^ diode 
transfers in "isolated" regime that leads to more complicated problem with free boundary. 

Thus we have the following main result of this paper. 

Theorem 4. Assume conditions (B 2 ) , (B 3 ), (B e ) and inequalities (14), (17), (23). Then 
the problem (I) possesses a positive solution P such that 

<Po > JxF((f , z 2 ), z 2 e [0, <p°] 

(<P°)"<j x F(<p° : z 2 ) : z 2 e[0 : <p°] ' 

% > G(j x ,z u a ), z x e [ip ,ip°] 
(a°)"<G(j x ,z 1 ,a°), z 1 eivo,¥°] ' 

where </?o — 5 2 x i ^ 3 is a lower solution of problem (Ai), ip° = a + (3x (a, (3 > 0) is an upper 
solution of problem (Ai) with condition Lpi > 5 2 ; a® = is a lower solution of problem (A4) 
with condition < a^ < J^>°(2 + <p°). 
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